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ABSTRACT 

SDSS 1257+5428 is a white dwarf in a close orbit with a companion that has been suggested to be a neu- 
tron star. If so, it hosts the closest known neutron star, and its existence implies a great abundance of similar 
systems and a rate of white-dwarf neutron-star mergers similar to that of the type la supernova rate. Here, 
we present high signal-to-noise spectra of SDSS 1257+5428, which confirm an independent finding that the 
system is in fact composed of two white dwarfs, one relatively cool and with low mass, and the other hotter 
and more massive. With this, the demographics and merger rate are no longer puzzling (various factors com- 
bine to lower the latter by more than two orders of magnitude). We show that the spectra are fit well with 
a combination of two hydrogen model atmospheres, as long as the lines of the higher-gravity component are 
broadened significantly relative to what is expected from just pressure broadening. Interpreting this additional 
broadening as due to rotation, the inferred spin period is short, about 1 minute. Similarly rapid rotation is only 
seen in accreting white dwarfs that are magnetic; empirically, it appears that in non-magnetized white dwarfs, 
accreted angular momentum is lost by nova explosions before it can be transferred to the white dwarf. This 
suggests that the massive white dwarf in SDSS 1257+5428 is magnetic as well, with B ~ 10^ G. Alternatively, 
the broadening seen in the spectral lines could be due to a stronger magnetic field, of ^ lO^G. The two models 
can be distinguished by further observations. 

Subject headings: white dwarfs: general — white dwarfs: individual (SDSS 1257+5428) 



1. SDSS 1257+5428 

An unexpected return from the Sloan Digital Sky Survey 
(SDSS; lYork et al. .2000) has been the discovery of white 
dw arf binaries wit h short periods, by Badenes et al. (2009j) 
and lMuUally etal.1 (120091) . The massively multiplexed SDSS 
spectroscopic observations consist of several 15 -minute in- 
tegrations^^ which are combined to yield the final spectra. 
iMullallv et al. (2009) and Badenes et al. (2009) took advan- 
tage of this approach and looked for rapid changes in radial 
velocity across the individual spectra. 

The first finds are already quite interesting. 
SDSSJ143633.29+501026.8 (orbital period, =1.15hr) 
and SDSSJ105353.89+520031.0 (Pb = 0.96hr) are both 
low mass (^0.3 M0) DA white dwarfs orbiting an unseen 
companion. For both systems, the measured mass function 
can be reasonably explained by invoking a white dwarf for 
the secondaryO The orbital periods are so short that the 
systems are expected to merge within a Hubble time. These 
two binaries increased the toll of such interesting double 
degenerates to seven (Mullallv et al. 2009). The number 
is likely to increase further; e.g., from a targeted se a rch of 
low-mass white dwarfs found in SDSS, iKilic et al.l (l2009h 
find four short-per iod binaries (including the two from 
iMuUally et al.l l20()9^. These authors also discuss in detail the 
fate of these systems: since the combined masses are likely 
below the Chandrasekhar mass, most will not explode as type 
la supernovae, but rather become AM CVn systems or R CrB 
stars. 

SDSS J125733.63+542850.5 (SDSS 1257+5428 hereafter) 

' On sabbatical leave from Department of Astronomy and Astrophysics, 
University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, 
Canada. 

^ We refer to the detected DA white dwaii' as the photometric primary (or 
"primary") and the unseen/fainter companion as the photometric secondary 
("secondary"). 



consists of a DA white dwarf primary in an orbit with period 
Pb = 4.6 hr and no measurable eccentricity. The inferred com- 
panion mass (as constrained by the spectroscopically inferred 
white-dwarf mass and the mass function) is above the Chan- 
drasekhar mass (in a probabilistic sense). If so, the companion 
is either a neutron star or a black hole, and, at the estimated 
distance, d w 5 pc, it would be the n earest ultra-compact ob- 
ject known; see lBadenes et al.l (120091) . 

The importance of the proximity of SDSS 1257+5428 be- 
comes apparent when we consider the distances to the near- 
est members of different sub-classes of neutron stars (see Ta- 
ble [TJ: the nearest millisecond binary pulsar (a binary con- 
sisting of a low-mass white dwarf and a pulsar and in a cir- 
cular- orbit), PSR J0437-4715; the thermally X-ray emitting 
middle aged neutron star, RX J 1856.5-3754; the nearest long 
period pulsar, PSR J2144-3933; the nearest ordinary pulsars, 
PSR JO 1 08-143 1 and PSR 0950+08 ; the nearest young pulsai-. 
Vela and the nearest 7-ray pulsar, Geminga. (For the ordinary 
pulsars, accounting for the usual beaming factor of 0. 1, results 
in a distance of about 130pc.) The distance c/ns ^ 150pc for 
the nearest neutron star(s) is in accord with the global demo- 
graphics of neutron stars as well as demographics of specific 
classes of neutron stars (based on birth rates). 

Specifically for SDSS 1257+5428, if its unseen companion 
is truly a neutron star then, given the short period, the sys- 
tem would arguably be a binary millisecond pulsar The lo- 
cal surface density of pulsars and millisecond pulsars is about 
30kpc"- (L vne et al.l il998). For millisecond pulsars the local 
density is between 30-45 kpc"^ (C ordes & Chernoft ,1997,) . 
The expected distance to the nearest millisecond pulsar is thus 
1 75-200 pc, consistent with the observations but inconsistent 
with the proximity of SDSS 1257+5428. 

Thus, a neutron-star companion to SDSS 1257+5428 would 
be very puzzling, unless one assumes that its proximity to the 
solar system is a statistical fluke. Motivated by this, as well 
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Table 1 






The Nearest Neutron Stars. 






Distance 




Name 


(pc) 


Ref. 



PSRJ0437-4715 

RXJ1856.5-3754 

PSRJ2144-3933 

PSRJ0108-1431 

PSR()950+08 

Vela Pulsai' 

Geminga 



156.3 ±1.3 

161!1« 

165! 

240+'f 

260+f 

286+|« 
250+^20 



Deller et al. 2008 

van Kerkwiik & Kaglm^OO? 

Deller et al. 2009 

Deller et al. 2009 

Brisken et al. 2002 

Dodson et al. 2003 

Fahertv et al. 2007 



as by t he alarming implications of a paper by lThompson et all 
(l2009h . we first reconsider the case for a neutron-star compan- 
ion. A reading of the Badenes et al. (2009) paper shows that 
it rests primarily on the high mass for the primary, with kine- 
matics offered as a supporting argument. We revisit these two 
issues in the next two sections, and then present observations 
showing that the system in fact hosts two white dwarfs. 

Before proceedi ng, we note that wh ile writing up our re- 
sults, a preprint by Marsh et al. (2010) appeared, which pre- 
sented lower-quality but much more numerous spectra, from 
which the presence of a second component is also evident. 
These authors discuss in detail the preceding and future evo- 
lution of this binary. In this paper, we focus on the demo- 
graphics, and on obtaining more accurate constraints on the 
white dwarfs in the system from our higher-quality spectra. 

2. REVISITING THE MASS OF THE COMPANION 

iBadenes et al.l (l2009t) present radial-velocity variations of 
SDSS 1257H-5428. The velocity data appear to be of high 
quality and so we accept the two inferences: orbital period, 
Ph = 4.555 hr and a velocity semi-major amplitude of the pri- 
mary, Ki = 322.7±6.3kms"'. The uncertainty in P is virtu- 
ally negligible. The mass function is thus 



M2 



Mj sin / 



PhKj 

(Mi+M2)2 2ttG 



: 0.66 ± 0.04 M, 







(1) 



where the uncertainty of 5.9% of A^2 arises from cubing Ki. 
The minimum mass of the unseen companion is AI2 (obtained 
by assuming that the system is seen edge-on and setting Mi = 
0)^ 

IBadenes et al.l (l2009h infer a mass for the primary based on 
fitting the Balmer absorption lines to model atmospheres of 
DA white dwarfs and find Mj = 0.92![J J^Mq. The large un- 
certainty in the inferred mass is due to the poor model fit to 
the data (more on this issue below). 

As noted by IBadenes et al.l ( 120091) . the unseen companion 
is quite massive (for a white dwarf) and exceeds the Chan- 
drasekhar limit for Mi > O.64M0. There is an 8% chance that 
Ml is less than this value and thus the above appears evidence 
for the conclusions of Badenes et al.. (2009). 

One issue is that empirically it has been found that spectro- 
scopic masses inferred for cool white dwarfs (Teff < 12, 000 K) 
are larger than the true masses. This has been a long-standing 
issue in the field of white dwarfs, but with no resolution in 
sight; see Tremblay et al. (2010) and Koesteretal. (2009a) 
for recent summaries. Spectral fitting yields P = g/ n, where 
P is the gas pressure, n is the mass opacity, and g is the local 
acceleration due to gravity. It has been suspected that in the 
cool white dwarfs, which have convective atmospheres, he- 
lium might be dredged up. This would pollute the atmosphere 



and decrease k, (relative to a pure hydrogen atmosphere). 
This hypothesis would then explain why the masses inferred 
spectroscopically are systematically higher for cool DA white 
dwarfs. Sensitive measurements in field white dwarfs, how- 
ever, have failed to detect the expected He I features at 5876 A 
tTremblay et al. 2010). 

Ignoring this very interesting and perplexing issue we pro- 
ceed in an empirical fashion. Using the run of inferred mass 
as a function of T^if for DA white dwarfs, we derive a mul- 
tiplicative correction factor of 0.9 for the mass inferred from 
spectroscopic modeling. The resulting decrease in mass to 
Ml = 0.8 ± 0.3 also decreases the mass estimate for the sec- 
ondary and thus weakens the case for a neutron-star compan- 
ion. 

Finally, and more importantly, the spectral fits to the data 
are poor This is acknowledged by Badenes et al. ( 2009), but 
it is worth stressing that such poor fits are not routinely seen 
in DA white dwarfs (e.g., see Tremblay et al. 2009). That typ- 
ical fits are good can also be seen from the small inferred un- 
certainties in Teff and \ogg. Taking the SPY survey of high 
sensitivity and high spectral resolution observations of white 
dwarfs as an example, and restricting ourselves to the sub- 
set of cool DA white dwarfs {T < 12, 000 K), we find that the 
largest uncertainty in \ogg is 0.03 (corresponding to 7% un- 
certainty) and the median 0.01 (corresponding to 2.3%). 

Thus, the poor fit for the Balmer Unes of SDSS 1257H-5428 
suggest that either SDSS 1257-1-5428 has a strange pathology 
(in which case modeling of the Balmer lines is fundamentally 
suspect) or that the companion is contributing to the light. If 
the latter, it might be responsible for some of the wiggles 
seen in the continuum, or, if it is a DA white dwarf, con- 
tribute to the observed B almer lines. F rom a more detailed 
look at the spectral fits of IBadenes et al.l (2009), one sees that 
these poorly match the narrow cores of lower Balmer lines 
and under-predict the higher Balmer lines. The presence of 
the higher Balmer lines is indicative of lower gravity, since 
these lines are very sensitive to pressure broadening (the size 
of a hydrogen atom is quadratically proportional to the upper 
level of the transition). Indeed, the fact that in the spectrum 
presented by Badenes et al. (2009), the line at HIO has depth 
comparable to that of H9 already informs us that \ogg < 7. 
This inference, combined with the wide wings shown by the 
lower Balmer lines, supports the idea that the observed spec- 
trum of SDSS 1257-1-5428 is produced by more than one DA 
white dwarf. 

3. REVISITING THE KINEMATICS 

Pulsars, both single and binary, are noted for their 
large space motion (ICordes & Chernof3 Il997h . The 
mean 3-D speed of the low-mass millisecond pulsars is 
84kms"'. IBadenes et al.l (120091) offer the kinematics of 
SDSS 1257-1-5428 as a consistency argument for a neutron star 
secondary. 

The proper motion of SDSS 1257H-5428 is 0.049masyr-'. 
At the stated distance of 50 pc the proper motion trans- 
lates to llkms"', which is not particularly distinguishable. 
Badenes et al. (2009), from their velocities, find that the mean 
radial velocity of the primary is 71 = -29 ± 5 km s~' . 

A compilation o f the mean radial veloc itiej^ of double de- 
generates listed by iNelemans et al.1 (12005!) shows a spread of 
-36 to 70kms~'. Correcting for the expected gravitational 

^ For consistency, we restricted the analysis to the radial velocity of the 
photometric primary. 
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Table 2 

Log of Observations and Radial Velocity Measurements 



Epoch 


Slit 


Cam 


r 


MJD 





V 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


12:13 


0.7 


B 


780 


0.51660 


0.86 




12:14 


0.7 


R 


340 


0.51470 


0.85 


-292 ±3 


12:21 


0.7 


R 


340 


0.51945 


0.87 


-262 ± 3 


12:28 


1.0 


B 


780 


0.52727 


0.92 




12:28 


1.0 


R 


340 


0.52491 


0.90 


-220 ±3 


12:35 


1.0 


R 


340 


0.52967 


0.93 


-173 ±3 


12:42 


1.0 


B 


780 


0.53709 


0.97 




12:43 


1.0 


R 


340 


0.53473 


0.96 


-120±3 


12:49 


1.0 


R 


340 


0.53950 


0.98 


-71 ±3 



Note. — Col. (1): Exposure start time on 2010 February 11. Col. (2): 
Slit width in arcseconds. Col. (3): Camera (B: Blue, R: Red). Col. (4): In- 
tegration time in seconds. Col. (5): Barycentric mid-exposure time relative 
to MJD 55238.0. Col. (6): Orbital phase calculated using the ephemeris of 
[Marsh et al. (2010), MJD 54845. 67470(8) + 0.18979154(9)£, where phase 
corresponds to inferior conjunction of the photometric primary. Col. (7): Ra- 
dial velocity calibrated relative to sky lines, and coirected to the solar system 
barycenter. There may be systematic effects larger than the formal uncertain- 
ties due to inaccuracies in our connection for the high-gravity component. 
For the observations, we used the Atmospheric Dispersion Corrector (ADC), 
and set up LRIS as follows. A dichroic at 5600 A separated the incoming 
light into a blue and a red channel. For the blue channel we employed the 
600 lines mm"' grism, blazed at 4000 A. The dispersion is 0.63 A per pixel 
and the spectral FWHM 4.1 A (for the 1 arcsecond slit). The detectors con- 
sisted of two 4096 X 2048 pixel Marconi CCDs with each 15 /im square pixel 
covering Of 135 on the side. For the red channel we used a grating with 
1200Unesmm"' and blazed at 7500 A. The dispersion was 0.4 A per pixel 
and the spectral FWHM 2.1 A (for the 1 arcsecond slit). The detector for 
this channel was a mosaic of two LBNL CCDs, of the same size and with 
the same plate scale as the blue ones (but read out binned by 2 in the spatial 
direction). 

redshiffl for the primary, using the masses as obtained from 
spectroscopic modeling, one finds that the true radial veloc- 
ities for these 24 double degenerates lie between -60 and 
50kms-i. 

For the corrected spectroscopic mass of 0.8 ± 0.3 M©, the 
gravitational redshift (converted to velocity) is 56;!;3Qkms"'. 
The l-a lower bound results in -55 ± 5kms"' - well within 
the expected kinematics of double degenerates. 

4. SPECTROSCOPIC OBSERVATIONS & DATA REDUCTION 

On the night of UT 2010 February 1 1, we undertook obser- 
vations of SDSS 1257H-5428 with the Low Resolution Imag- 
ing Spectrograph (LRIS) mounted on the Keck I 10-m tele- 
scope (Oke et al. 1995). The night was clear and the seeing 
typically about 0."7 (red camera) and 0."8 (blue camera). The 
observations are summarized in Table |2] We started a se- 
quence of observations of SDSS 1257H-5428 using a 0."7 slit 
and later changed to 1"0 to better match the seeing. Fol- 
lowing this, we observed a number of internal light sources 
(incandescent for flat fields, Hg/Kr/Zn/Ar and Ne/Ar lamps), 
as well as, at a later time, observations of the flux standard 
Feige 34 (Oke 1990) with the same spectrograph setup. 

For the reduction, we used the European Southern Observa- 
tory (ESO) Munich Image Data Analysis System (MIDAS), 
and routines running in the MIDAS environment. For all im- 
ages, we subtracted bias as determined from the overscan re- 

* Throu ghout this article we use the mass-radius fitting formula of 
INauenbe r? 1 1 972) which provides a good approximation to the classical re- 
sults of Hamada & Saloetei ( 1961). The gravitational or Einstein redshift, 
converted to velocity, is = 0.633M/Skm s"' where M is the mass and R 
is the radius. 



gions. For the blue images, we subsequently corrected for 
small-scale variations in efficiency by dividing by a spatially 
averaged flat field, normalized using a third-degree polyno- 
mial, and with the bluest, poorly exposed part short-ward of 
4000 A replaced by unity. For the red images, we simply di- 
vided by the flat field, normalized using a bi-linear fit. We 
extracted the spe ctra using an optimal extraction procedure 
similar to that of iHornd (Il986 i). after subtracting sky deter- 
mined from neighboring regions. 

Wavelength calibration was done using arc spectra. For the 
blue arm, fifth degree polynomial fits were required to give 
adequate dispersion solutions, with typical root-mean-square 
residuals of 0.02 A (for 14 lines). For the red arm, a fifth- 
degree polynomial gave residuals of ^ 0.017 A (for 27 lines). 
We used the same wavelength solution for all spectra, but cor- 
rected for possible drifts using the oxygen sky emission lines 
at 5577.340 and 6300.304 A, in the blue and red, respectively. 
For the blue, the line is at the edge of our wavelength range, 
and hence our wavelength scale may be slightly off. Indeed, 
we find that the velocities inferred from the blue spectra show 
a slit-size dependent offset from those inferred from the red 
spectra; hence, we will use the blue side only to (attempt to) 
measure relative velocities. 

For flux calibration, we first corrected all spectra approxi- 
mately for atmospheric extinction using a curve made b y com- 
bining the CFHT values (B eland. Boulade. & Davidgen 19881) 
short-ward of 5200 A with the better sampled La Silla val- 
ues long-ward of 5200 A (ESO users manual 1993; see also 
'Tug!ll977'). Next, for the blue spectra, we calculated response 
curves by comparing ou r obs e rved s pectra for Feige 34 with 
the calibrated specti-a of lOkd (Il99d) as provided by STSclfl 
we slightly smoothed our spectra to match the resolution, di- 
vided the two, and smoothly interpolated the ratio. Since 
all spectra were taken through a relatively narrow slit, this 
should give good relative calibration, but the absolute calibra- 
tion may be off. Folding our flux calibrated spectra through 
the SDSS M, g, and r response curvesQ we infer g = 16.6, 
M-^ = 0.41, and g-r = 0.1 1, which i s in go od agreement with 
the values listed by iBadenes etaTI ( 120091) : 16.844 ± 0.004, 
0.511 ±0.010, and 0.112 ±0.007, respectively, especially 
considering that our spectra do not completely cover the u 
band. For plotting purposes, we rescaled our blue and red 
spectra to match the g and r, magnitudes, respectively. 

5. SPECTRAL ANALYSIS 

The spectra show Balmer lines up to H12 and distinct asym- 
metries in some of the Balmer absorption features, with the 
red wing shallower than the blue one (Figure[Tl note that there 
is no sig n of wiggles in the c ontinuum such as those seen in 
Fig. 5 of'B adenes et al.ll2009l) . This suggests that two sources 
contribute to the Balmer features, with the second red-shifted 
relative to the one responsible for the line cores. The asym- 
metry is most pronounced in HS and He, and less so for both 
lower and higher Balmer lines. The former suggests that the 
temperature of the second component is higher than that of 
the most prominent one (but that the emitting area is smaller), 
so that it contributes mostly at bluer wavelengths, while the 
latter suggests its gravity must be higher, so that the higher 
Balmer lines are less strong. 

Granted that there are two sources contributing to the 

^ ftp://ftp.stsci.edu/cdbs/calspec/ 

* http://www.sdss.0rg/dr6/instruments/imager/index.html#filters 
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Figure 1. The Balmer Series of SDSS 1257+5428. The spectra are ordered chronologically from bottom to top, and offset by 0.2 units (0.14 units for Ha). The 
second trace has offset 0. The presence of Balmer lines up to H12 shows that one component of the binary has low gravity, while the wide wings on the lower 
Balmer lines suggest a high-gravity companion, with the asymmetry relative to the cores due to orbital motion and differences in gravitational redshift. The 
lower, red curve shows the model fit described in the text, and shown in more detail in Figure[2] 



Balmer series we were, nonetheless, puzzled by the absence 
of double line cores. After some experimentation and contem- 
plation we invoked rotational broadening for the secondary 
star. 

5.1. Model Fining 

The above qualitative impressions are born out by fits with 
model atmospheres. For these fits, we modeled the blue 
spectra, in the range of 3700-5200A (i.e., the range show- 
ing lines), with a combination of two hydrogen model atmo- 
spheres, taken from a set kindly provided by D. KoesterQ 

We scanned a grid in effective temperature and surface 
gravity for each component, as well as a set of rotational ve- 
locities < VrSin; < 2000km s"' for the higher gravity com- 
ponent (in steps of 200kms ~'; we conv olved the model us- 
ing the analytical profile of iGravl (12005) . with a limb dark- 
ening coefficient of 0.^|). For the lower-gravity component, 
we took the smearing due to chan ging orbital velocity into 
account (using the orbit of lMarsh"e t al. 2010; the largest ef- 
fect is 80kms"', substantially smaller than our resolution; the 
motion of the high-gravity component during an integration is 
negligible). 

We accounted for the spectral resolution by convolving all 
models with a truncated Gaussian (with width and truncation 
matching the seeing and slit width, respectively), and resam- 
pling them at the observed wavelengths. At each grid point, 
we fitted for the two best-fit velocities, the flux ratio, and a 
normalization (which, to account for possible errors in our 
flux calibration, we allowed to depend quadratically on wave- 
length). 

Our fits yield a reasonably well-defined minimum, at the 



Table 3 

Spectroscopic Properties of the White Dwaii's in SDSS 1257+5428 



Property 


Phot. Primary 


Phot. Secondary 


Fit parameters 






Teff (K) 


6250 ±250 


13000 ±500 


logg (cms"2) 


6.0 ±0.3 


8.5 ±0.2 


V,- sinr (kms"') 


<100 


800 ±300 






0.21 ±0.03 


Derived parameters 






M(Mq) 


0.15 ±0.05 


0.92±0.13 




0.042 ± 0.008 


0.0089 ±0.0012 



parameters listed in Table 
whole spectrum in Fig. |2 



3 We show the fit overlaid on the 
One sees that, overall, the fit is 



' This is Koester's most recent version of grid of white dwarf spectra. 
The basic model is described in Koester et al. 1 2009b) and the latest version 
includes improved treatment of pressure broadening (Tremblay & Bergerq^ 
l2009h . In detail, the giid spans 6000 < Teff < 30000 K in steps of 250 
to 1000 K at low and high temperature, and 6.0 < logg < 9.5 in steps of 
0.25 dex. 

* The limb darkening coefficient is inferred from mo del specific intensitie s 
used previously to fit spectra of pulsating white dwarfs ^ Clemens et alJ2000ll . 
We ignored the variation in limb darkening over the lines. 



Note. — The parameters are for a model in which the binary is composed 
of two white dwaii's, with the secondary rotating rapidly. The limit to the 
rotation of the primary is derived separately from Ha; at the limit, the fit to 
the core is obviously worse. Note that the model fit is not formally accept- 
able (see text); the quoted uncertainties are what we believe are conservative 
estimates. 



good, also outside the fitted range (indeed, it reproduces the 
spectral slope remarkably well). 

In detail, however, the fit has problems: the highest Balmer 
lines are predicted slightly too weak, and some line cores are 
not reproduced accurately (Ha being not deep enough, and 
H7 and H8 too deep). These effects might indicate that the 
true gravity of the primary is slightly lower than that in our fit 
(which is at the lower boundary of our grid). 

Furthermore, the red wings of some of the lines are poorly 
matched: that at H7 is fainter than observed, and those at 
He and H8 are too bright. We believe this difference is real 
rather than due to, e. g., inac curate flux calibration. Indeed, 
in the spectra of Ma rsh et al.l (12010), the same discrepancies 
are se en. One possible reason, which we discuss further be- 
low 06. 5b . is that the higher-gravity white dwarf is magnetic 
and that its line profiles are (partially) broadened by Zeeman 
splitting. 

Because of the above problems, even though the fit is 
hugely supe rior to one using on ly a single white dwarf (cf. 
Figure 5 of iBadenes et al.l [20091) . it is still not formally ac- 
ceptable; for the spectra taken through the 1" slit, we find 
1 .7 (for 2408 data points and 13 parameters; Xred —1-5 
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Figure!. Spectrum of SDSS J1257+5428 and a model decomposition in two wliite-dwarf spectra. In tlie left-liand panel, the top curve shows the observed 
spectrum, offset by 0.4 units (with the location of the dichroic indicated; we used the second set from Table[2] with two red spectra averaged). The spectrum is 
repeated below, with the model overlaid (fit to the 3700-5200 A range [blue points]). The two red lower curves show the two components in the model, one that 
has low gravity and relatively cool temperature, and one that has higher gravity and temperature, and is rotationally broadened. In the right-hand panel, details 
of the fit around the Balmer lines are shown, with wavelength shifted to the rest frame of the low-gravity component, and offsets of 0. 1 in flux added between all 
profiles beyond H/3. 



for the single, less well-exposed spectrum taken through a 0"7 
slit). While the fit is internally consistent, in that the three 
blue spectra give the same best-fit grid point, and that, e.g., 
fits to Ha and to higher Balmer lines give roughly consistent 
results (for the latter, though, there is a much larger covariance 
with the properties of the high-gravity component), we cannot 
know what offsets are caused by the systematic differences. 
We give what we believe are conservative error estimates in 
Table E] 

Finally, we searched for the presence of He I A5876 A (see 
Inspecting the spectrum, we determined that a 2% de- 
pression would be clearly noticeable. If an absorption feature 
existed at the limit of instrumental resolution, then the equiv- 
alent width limit would be 0.02 x 3 A=0.06 A. If the feature 
was broadened by SOOkms"', corresponding to AA ^ 15 A, 
the equivalent width limit is 0.3 A. 

5.2. Comparison 

Comparing our inferred parameters with those given by 
IMarsh et all dloTol) . one finds qualitative agreement, in in- 
ferring a cool, low-gravity and a hotter, high-gravity com- 
ponent, but what would appear significantly different param- 
eters (MarsJietalJ 12010 find r^ffi = 7200±350K, loggi = 
6.85 ±0.10, reff,2 9800± lOOOK, and log^2 ^ 9.0 ±0.4). 
It appears that the most important difference is that in Teff a: 
if we fix this temperatu re, the remaining p arameters become 
comparable to those of IMarsh et al.l (120101) . though the fit is 
much worse (x^d — ^■^)- Ths same does not hold if we 
fix other parameters (e.g., for Teff.i = 7250, the fit quality is 
similarly poor for a large range in Teff a)- The underlying 
problem may be that the Balmer lines have similar strength 
at 10,000K and 13,000K. This leads to a de generacy for 
normalized spectra such as those presented bv Marsh et al.l 
( 12010 ). Usually, such degeneracies can be broken by opti- 
cal photometry, but here one can change the temperature of 
the other component t o compensate. Fort unately, the ultravi- 
olet data presented by IMarsh et"an (120101) show conclusively 
that the high-temperature solution is the correct one. 

Turning now to the velocities, Uke Marsh et al.i (120101) . we 



find that our fits yield precise measurements of the motion of 
the photometric primary, in both our blue and red spectra, but 
not of that of the secondary. For the primary, for the blue 
spectra the accuracy is limited by flexure in the spectrograph, 
which we appear to be able to correct only partially using the 
5577 A sky line. In Table |2] we list the barycentric veloc- 
ities inferred from Ha; these are consistent with the radial- 
velocity amplitude of 330kms~' and the systemic velo city of 
28.9 k m s~' found by iBadenes et aH (12009) and .Marsh et al.1 
(120101) . For the secondary, we find that while for every choice 
of parameters, its velocity can be measured in the blue, and 
is consistently found at a positive offset from the primary (as 
expected given the orbital phases and its larger gravitational 
redshift), the precise values depend strongly on the tempera- 
ture of the primary, with changes of 250 K leading to velocity 
differences of ~'50, 40, and 30km s"' for the three different 
spectra. Thus, we cannot measure its motion. 

5.3. Inference 

From the best-fit model parameters for the second ary (Ta- 
ble |3]l and using the cooling models of ' WoodI (Il992h . we de- 
duce the following: M2 = 0.92 ± O.OMq and R2 = 0.0089 ± 
0.0012/^0, and a cooling age of t, = 0.7 ± 0.3 Gyr. 

For the primary, we start with the normalization in the spec- 
tral fits, which yields the ratio of emitting areas for the two 
components, and thus R2/R1 = 0.21 ±0.03. Combined with 
the estimate for R2, one finds Ri = 0.042 ± 0.008/?©. The 
spectral fits yield loggi = 6.0 ±0.3 (but more about this be- 
low). The mass of the primary can be directly computed from 
loggi and the radius, Ri. The mean value is 0.07Mo; how- 
ever, due to the large error in \ogg, the distribution of Mi has 
a long tail with a variance comparable to that of the mean. 

Thus, we cannot determine the mass directly, but we can 
compare the parameters with those for other low-mass white 
dwarfs. The first w ere seen as companions t o millisecond 
pulsars (for a review, Ivan Kerkwiik et alll2005l) and more re- 
cently field objects were found via SDSS dKilic et alj [2007). 
Some notable objects are SDSS J0917H-46 {logg = 5.48 ± 
0.03, r.ff = 1 1288 ± 72 K; .KiUc et al.. .2007.): NLTT 11748 
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(logg = 6.20 ± 0.05, Teff = 8540 ± 50K; iKawka & Vennesl 
20^ ; LP 400-22 (logg = 6.35 ± 0.05, Teff = 1 1 170 ± 90 K; 



Vennes et alj|2009l) . The latter is also notable because it has 



a high space motion, even though its companion is another 
white dwarf. 

From evolutionary considerations, the true mass of the pri- 
mary in SDSS 1257H-5428 is likely higher than the above es- 
timate of 0.07 Mq from logg, since given its current orbital 
period, the progenitor of must have evolved off the main se- 
quence, which implies it must have had a core mass of at least 
^O.15M0. This would imply a true log^i > 6.3. This is 
on the upper limit of our measurement, but we note that the 
model atmospheres appear not to entirely reliable in this low 
gravity regime. Specifically, models from different groups 
gave results different by 0.4 dex for the l ow-gravity, logg ^ 
6.5 c ompanion to PSR J1012H-5307 (see Ivan Kerkwi^k et alj 
l2005h . 

6. SDSS 1257+5428, A DOUBLE DEGENERATE: RAMIFICATIONS 

Above, we estabhshed that SDSS 1257H-5428 is not a neu- 
tron star binary but a bona fide double degenerate system with 
parameters given in Table |3] Here, we explore the conse- 
quences of this conclusion. 

6.1. Distance & Demography 

The distance to SDSS 1257H-5428 is increased because the 
flux of SDSS 1257H-5428 is now dominated by a lower mass 
white dwarf (with a larger radius and hence a larger an- 
gular diameter relative to that in the iBadenes et al.l 120091 
model). From our fit (see Table O, we infer a distance to 
SDSS 1257H-5428 of about 140 pc (instead of 50pc). 

Is SDSS 1257H-5428 unique? There exists 129 systems (sin- 
gle white dwarfs, white dwar f binaries and do u ble degenerate 
binary systems) within 20 pc jSion et al.l2009l) . iHolberg et alj 
( 120081) find seven double degenerate systems within 20 pc. 
Thus, within 140pc one expects to find about 2400 double 
degenerates. A different estimate can be made by combin- 
ing the local d ensity of white dwarfs o f 4.5 x lO'-'pc"^^ with 
the finding of iMaxted & MarshI (1 19991) fliat about 1.7-19% 
(95% confidence level) of white dwarfs are double degener- 
ates. This yields an expected 900-10000 double degenerates 
within 140pc. Thus, at first blush, SDSS 1257H-5428 is not a 
rare system. 

Furthermore, the systemic radial velocity of 
SDSS 1257H-5428 is now 71 - 7^(1), where the former 
is the mean radial velocity of the primary (see and 7£(1) 
is the gravitational redshift for the primary. With the revised 
mass, the redshift is small, 7^(1) =2 to 5kms"' (for Mi =0.1 
to 0.2Mq). Thus, the 3-D velocity of SDSS 1257H-5428 is 
quite consistent with that of other double degenerates. 

The (exponential) scale height for the double degenerates, 
given their measured velocity dispersion, is similar to that of 
(say) G dwarfs, or about 350 pc. 

6.2. Coalescence Time-scale 

Double degenerates, especially in compact binaries, lose 
energy by gravitational radiation and eventually coalesce. The 
time for the orbital period to halve is (iLang,1980„ p. 600): 



]_dP^ 
Pb dt 



96 G AiM^ 

1 

2.8 X lO^yr 
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0.3 0.4 

cos(?) 

Figure 3. Mass of the photometric secondary as a function of cos i where 
/ is the angle between the orbital angular momentum and the line of sight 
(from bottom to top, Mi = 0.10. 0.15, 0.20. 0.25 M(7)). Here we ad opt the 
radial-velocity ampHtude for the primary from Marsh etaTI I2010D , Ki = 
330 ± 2kms"' . The pair of dotted lines represent the uncertainty in M2 (for 
a given value of Mi ) arising solely from the 1.8% uncertainty in the mass 
function, M2- The nominal value of Mi = 0.15Mq is shown by the thick 
line. The bounds on the secondary mass M2 from spectral fitting is shown by 
the horizontal dotted lines. The totality of the data argue that M| < O.2M0 . 

where M = Mi +M2 and = M1M2/M. Integrating this equa- 
tion to zero period yields a coalescence time-scale of 



TGw= 1-05 X 10^ yr 



M 



-2/3 



Mq 



Ihr / 
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(3) 



For Ml = 0.15 ± 0.05 M© and M2 = 0.92 ± 0.13 Mg we find 
"^Gw = 4.4Gyr. At the 1 percentile level th e range is 2.5- 
15Gyr. Incl udin g the cooling time-scale (iPhinneyl 119911) 
noted above 05.3l l results in a total duration of ~5Gyr. 

6.3. Merger Rate 

[Thompson et alj (120091) . ac cepting the masses o f and dis- 
tance to SDSS 1257H-5428 of I Badenes etal] (12009 1). found a 
very large rate of coalescence of SDSS 1257-i-5428-like sys- 
tems - comparable to that of la supernovae. The revision of 
the distance, masses and scale height reduces the coalescence 
rate by a larg e fac tor: using a scale height of 0.5 kpc instead 
of 4kpc (see ^61] ); a distance of 140 pc to SDSS 1257H-5428 
instead of 5 pc 06.11 ): and a lifetime of ^5Gyr instead of 
2Gyr 06.21 ). our estimate of the Galactic coalescence rate is 
a factor 4/0.5 x (140/50)^ x 5/2 smaller, or 



7mw ' 



1 X 10"''A^yr" 



(4) 



where = 1 is the number of such systems detected (and con- 
sequently suffers from severe Poisson uncertainty). 

The rate given in Equation H] is not alarming (correspond- 
ing to a mean time between events of ^ 1 Myr). Furthermore, 
given that Mi -I-M2 does not incontrovertibly exceed the Chan- 
drasekhar limit, the outcome of this coalescence is (conserva- 
tively) not a type la supernova explosion, but rather the for- 
mation an R CrB heli um giant or an AM CVn system (see the 
discussion in lKilic et al. 2009). The rate given in Equation|4] 
is small enough that, even if the combined mass exceeds the 
Chandrasekhar limit, SDSS 1257-i-5428-like systems are not 
major contributors to the supernovae la rate. 

6.4. Broadening: Rotation ? 
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As noted in ^ we are forced to invoke broadening for the 
Balmer series of the secondary, and modeled this as rotation. 
We consider this simplest possibility first. Our inferred rota- 
tion rate is v,.sin/ = 800± 300kms~' (see TableO, where / is 
the inclination between the rotation axis of the secondary and 
the line of sight. We will assume that the spin of the secondary 
is aligned with the angular momentum of the binary. 

For the entire range of plausible values for M2 and the 
bounds for Mi we find that cos / < 0.4 and thus sin/ > 0.84 
(see Figure|3]l. In view of the errors on vsin/, we will approx- 
imate sin / « 1 . Using our derived value for the radius of the 
secondary, R2 = 0.0089 ±O.OO12/?0 (see ®, we find that the 
rotational period, P ~ SO^jg s, where the 30% uncertainty in 
vsin / dominates the uncertainty. We will adopt a round fig- 
ure of 60 s for the rotation period of the secondary, and write 
corresponding angular frequency as uj = 2tt/ P. 

6.4.1. Spin-up by Accretion 

The simplest possibility that comes to our mind as an ex- 
planation for rapid rotation of the secondary is spin-up by ac- 
cretion of matter In order to lead to significant spin-up, this 
should happen in a relatively long-lived phase, and hence it 
seems unlikely the system was brought to its present state by a 
common-envelope phase (furthermore, to survive a common- 
envelope phase, the system would have to have been rela- 
tively wide, and one would expect a somewhat more massive 
helium-core white dwarf). Instead, it seems most likely the 
evolution would be similar to what is invoked to explain the 
millisecond periods of neutro n stars in binaries for both long 
and short periods (Ph > 1 d, IWebbink et all [19831: Ph < Id, 
iPylyser & Savonijgll988h . Depending on the evolutionary 
status of the mass donor and the orbital separation, mass trans- 
fer to a higher mass accretor (in this case the white dwarf 
photometric secondary) is driven by nuclear evolution, loss 
of angular momentum via stellar winds ("magnetic braking") 
or gravitational wave radiation. The typical accretion rate for 
short-period binaries is ^ 10~^Moyr"', corresponding to an 
accretion luminosity GM2M/R2 ~ 10^"* ergs"'. 

We first consider the case where the accreting white dwarf 
is unmagnetized. Its moment of inertia is I2 = ^2^2^?, with 
k2 = . 143 appropriate for a « = 1 .5 polytrope dBrooker & Olid 
[19551) . For an unmagnetized white dwarf, the accreted matter 
has the specific orbital angular momentum appropriate for the 
radius of the white dwarf, or ^2^2^, where the orbital angular 
frequency 0.2 = {GM2/R2^y'^ ~ 7s (for M2 w IMq). Thus, 
for an accreted mass AM, we have 



/\Mn2R2 =k2M2R2 
AM UJ 
=k2—=Qmi. 

M2 ^2 



(5) 



It is empirically known that accre ting white dwarfs in CVs 
do not rotate rapidly (e.g.. lGansicke et al. 2005). and one de- 
duces that the accreted angular momentum is lost, presumably 
in nova explosions. This requires that for most CVs there is 
poor angular mo mentum coupling betw een the accreted mat- 
ter and the core (Livio & Pringle 1998). Each nova cycle, in 
effect, leaves both the mass and the spin state of the accreting 
white dwarf unaffected. 

For the low accretion rates one expects for the progeni- 
tor of SDSS 1257H-5428, nova explosions will occur These 
happen when the accreted mass exceeds a certain "ignition" 
value, the value of which depends on the mass of the white 



dwarf. According to iTownslev & Bildstenl (120041) the igni- 
tion masses are Mig « (6, 3, 2) x IO'^Mq for (0.6, 1 .0, 1 .2)M0 
white dwarfs. Thus, assuming that the white dwarf accumu- 
lated no more than Mig following the very last nova explosion, 
the spin period can be no faster than about an hour 

The minimum amount of accreted matter can be reduced 
for the case when the accreting white dwarf is magnetized. 
In this case, the accreting matter has the specific angular mo- 
mentum \/GM7a, where oi at the Alfven radius. This matter 
has more specific angular momentum compared to the previ- 
ous case by the factor y^T^T^, which, once the white dwarf is 
spun up such that rA approaches the corotation radius, scales 
as cx (fl2/ujy^^ via Kepler's third law. Hence, in this case the 
minimum accreted matter is reduced to 



AM 



■■k2 



02 



1/3 



= 0.008M, 



0- 



(6) 



However, even this mass exceeds the ignition masses by more 
than an order of magnitude. Thus, we reject the spin-up-by- 
accretion model, at least for the case where the white dwarf 
is not magnetized or magnetized sufficiently weakly that the 
accreted envelope does not couple to the core0 

6.4.2. Intermediate Polar Model 

There exists a class of cataclysmic variables called interme- 
diate polars whose white dwarfs are found to be spinning quite 
rapidly. These consist of a strongly magnetized white dwarf 
accreting matter from the companion. The most famous ex- 
ample is AE Aquarii with a spin period of 32 s. A number 
of intermediate polars have spin periods of about a minute 
Motivated by the similarity of the spin periods of intermediate 
polars and that inferred for SDSS 1257H-5428, we explore an 
intermediate polar model for SDSS 1257H-5428. 

Intermediate Polars undergo nova explosions (e.g., GK Per 
= Nova Persei 1901 ; DQ Her = Nova HercuHs 1934; see Livi3 
FT983I) . (Novae also have been seen from polars, in which 
the white dwarf rotation is magnetic ally locked to the orbit; 
e.g. y i500 Cyg = Nova Cygni 1975; IChlebowski & Kaluznvl 
Il988t ) Thus, the rapid rotation seen in intermediate polars, 
despite the novae, shows that the white dwarfs manage to re- 
tain the accreted angular momentum. We infer that the strong 
magnetic field ensures rapid coupling between the accreted 
matter and the rest of the white dwarf. Subsequent novae ex- 
plosions merely result in ejection of the matter and a very 
modest fractional loss in angular momentum, |Mig/(A:2M2) 
(where we assumed /gj = ^Mi„R\ for the ejected shell). 

Following the accretion of the minimum mass, the white 
dwarf will not be spun up any further, but rea ch an equilibrium 
spin period, given bv (iGhosh & Lambl[T979l) . 



P,^ = 6.3 X Bfl'l^^ 



M 



R 



5 X 10*^ cm 



(7) 



Here, B = 10^B4G is the dipole field of the white dwarf and 
L = 10^'*L34ergs~' is the accretion luminosity during the ac- 
cretion phase. With M = O.9M0 and R = 0.009 Rq ~ 6 x 

' It may be possible to find evolutionary histories in which mass transfer 
is much more rapid, or that have, e.g., an initial phase with rapid mass trans- 
fer, when the donor is still more massive than the white dwarf. We did not 
investigate this further 

See http://asd.gsfc.nas a.gov/Koji.Mukai/iphome/iphome.htmll for a cat- 
alog of Intermediate Polars. 
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10** cm, we see that a modest (dipolar) field of 1.3 x lO^G 
is sufficient to account for the inferred period of ^60s. 
The dynamical equation for the spin-up is given by 

^ ~ 3 X lO-'^Hzs-' L"\'"r^ ( ^ V^' 

2^^3x10 Hzs ^34M30^5o(^^^j (^5^^08emj ' 

(8) 

where /i = BR^ is the magnetic moment and the normalization 
is that for a 10'* G dipole field and radius appropriate for a 
solar mass white dwarf. The white dwarf would be spun up 
(assuming a suitable accretion rate) in less than two million 
years. 

We conclude that the intermediate polar model is a good 
explanation for the origin of spin in SDSS 1257+5428. As 
an aside we make the following observation. The fastest in- 
termediate polars have periods of about a minute. This then 
empirically suggests that a magnetic field strength of ^ lO^G 
is the minimum field strength required to rapidly couple the 
white dwarf and the accreted envelope. 

6.5. Broadening: Magnetism? 

Above we assumed that the observed broadening in the 
high-gravity spectrum is due to rotation. An entirely differ- 
ent possibility is that it is magnetic. In a strong magnetic 
field, the Balmer series are split into three components: two 
senses of circular polarization and an unpolarized component. 
To first order, the circular polarization components are shifted 
by half the electron cyclotron energy, or 0.0057 eV (B/IO^G), 
on each side of the unperturbed unpolarized component. At 
second order, all lines will shift bluewards, with the shift sub- 
stantially stronger for transitions to higher excited states. 

For the field strength of lO^G inferred above, the total ex- 
pected spread is 0.001 1 eV, which corresponds to a "velocity 
width" of about 200 km s"' at Ha (and smaller proportionally 
with for the higher Balmer lines). This in itself will not 
have a noticeable effect on the model fits. 

However, a possibility is that the white dwarf is not ro- 
tating rapidly at all, but has a stronger field strength, say 
approaching 10^ G. In this case, the observed broadening 
could be almost enti rely due to Zeem an splitting (e.g., as in 
WD 0637+477; Sch midt et al.llT992l) . Indeed, this might be 
the cause for the relatively poor spectral fits. We looked for 
the expected decrease in broadening and increase in blueshift 
for the higher Balmer lines, but we were unable to conclu- 
sively accept or reject the hypothesis that the Balmer lines are 
broadened by Zeeman splitting. 

7. TESTING THE MODEL 

We found that SDSS 1257+5428 is composed of two white 
dwarfs, one relatively cool one with low-gravity, and a higher- 
gravity, hotter one, whose lines are substantially broadened. 
A question that is left is what is the cause for the broadening, 
rapid rotation with a period of about a minute, or a strong, 
^ lO^G, strength field. The clearest way to distinguish the 
two models would be by spectro-polarimetry. For the case 
of rapid rotation, this would also allow one to test whether, 
as we contended, such rapid rotation requires the presence 
of a weaker field, of ^ lO'^G, similar to the fields inferred 
for intermediate polars. Rapid photometry (and spectrometry) 
could perhaps even reveal the rotation period. With \ogg = 9.5 
and Teff = 1 3 , 000 K, the secondary whit e dwarf lies close to 
the ZZ Ceti strip dGianninas et al.ll2006[) and so a search for 
pulsations could be quite productive. 



Apart from its intrinsic interest, SDSS 1257+5428 may also 
become a useful test case for white dwarf models. With suf- 
ficient sensitivity and orbital coverage, it should be possible 
to measure velocity curves for both components (i.e., derive 
M2/M1) and determine the difference in gravitational red- 
shift, (i.e., (M/R)2-(M/R)i). Further model atmosphere fits 
would yield much more accurate temperatures and gravities 
(i.e., (M/7?^)i.2), and a precise ratio of the emitting areas (i.e., 
{R2/R\)^)- Combined with the theoretical mass-radius rela- 
tion, the system is thus strongly overconstrained, and can be 
used to test various assumptions. This would be espec ially 
valuable for the low-mass companion, since, as noted in 35.31 
for the very low-mass white dwarfs the atmosphere and mass- 
radius relation are currently not as securely established as is 
the case for regular white dwarfs. 

We thank E. S. Phinney, G. Nelemans, C. Badenes, and L. 
Bildsten for discussion. LRIS has been upgraded since its 
commissioning and as a result is perhaps now the most effi- 
cient single object optical spectrometer. We thank the teams 
which made these improvements possible (J. McCarthy & C. 
Steidel, leaders for LRIS-Blue upgrade; C. Rokosi, leader for 
the LRIS-Red upgrade; and J. Miller & D. Phillips, leaders 
for the ADC sub-system). We are grateful to the staff of the 
WM Keck Observatory for their excellent service and to the 
librarians who maintain the ADS and Simbad databases. 
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